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Abstract. The cosmic-ray propagation code GAL- 
PROP has recently been enhanced in the areas of 
Galactic magnetic field and synchrotron radiation. 
A full 3D treatment is now implemented. We de- 
scribe these enhancements and show applications 
to the interpretation of non-thermal radio surveys 
and gamma rays, to probe the Galactic cosmic-ray 
electron spectrum and magnetic fields. 
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I. Introduction 

Cosmic-ray electrons and positrons produce syn- 
chrotron radiation by gyrating in the magnetic field of 
our Galaxy. The same electrons are also partly responsi- 
ble for the Galactic diffuse gamma-ray emission, in the 
hard X-ray to gamma-ray energy range Q]. Since cosmic 
rays, gamma rays, synchrotron emission and Galactic 
magnetic field are interrelated, combining such observa- 
tions constrains cosmic-ray electrons and components of 
the interstellar medium (ISM) responsible for the diffuse 
emission over this broad energy range. 

At the conference we will present a study of the 
synchrotron emission in order to constrain magnetic field 
models of the Galaxy, using the knowledge from the 
gamma-ray and electron data of the Fermi telescope. 

We use the GALPROP code for the calculation of 
Galactic cosmic-ray propagation (O, Q, Q], 0, (6) 
and (3). The code is developed to reproduce direct 
cosmic-ray measurements, gamma rays and now syn- 
chrotron radiation data. More details of the GALPROP 
code and its developments are given in [8] and in the 
next paragraph. 

II. GALPROP DEVELOPMENTS 

Synchrotron emission depends on the model of the 
magnetic field. Hence, an implementation of a proper 
3D model of the Galactic magnetic field is fundamental 
to have a good description of the diffuse synchrotron 
radiation and its latitude and longitude distribution. Until 
now, a simplified magnetic field model in 2D was 



implemented in GALPROP, using only a random com- 
ponent with an exponential decrease in Galactocentric 
radius and height above the plane characterised by a 
scale length in each dimension. With the availability of 
excellent radio continuum surveys from tens of MHz to 
tens of GHz, including the WMAP satellite data, more 
sophistication is desirable. Hence GALPROP has been 
extended to include a general model of the B-field and 
the associated synchrotron emission. We have introduced 
full 3D models for both regular and random magnetic 
fields. The routines have been adapted to compute the 
emission in 3D; with the previous 2D formalism the 
asymmetries in the synchrotron emission were not de- 
scribed. 

At present the Galactic magnetic field is not well 
known. Hence many models with different configura- 
tions can be found in the literature. For our study we 
chose some of them as a starting point and made our 
own choice of the parameters. The formulation of the 
regular component of the magnetic field introduced in 
GALPROP follows a spiral form. The regular component 
is modelled as a vector field B reg , the random field 
as a scalar field B ran (x,y,z), with the option also to 
represent it as a vector in case explicit modelling of 
random orientations is required. 

Skymaps, longitude and latitude profiles at a variety 
of energies/frequencies, as well as spectral index plots 
are produced as outputs so that the diffuse emission for 
different models of the magnetic field can be compared 
with data. 

With the present version of GALPROP only the 
total synchrotron intensity can be compared with data; 
polarized emission will be implemented in future work. 

Electrons and positrons lose energy by inverse Comp- 
ton and synchrotron radiation, and the latter is included 
in GALPROP self-consistently using the total field of 
the adopted model. For the spectrum of particles (here 
electrons and positrons) computed by GALPROP at all 
points on the 3D grid, we integrate over particle energy 
and line-of-sight to get the synchrotron intensity for the 
regular and random fields. The resulting synchrotron 
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skymaps for a user-defined grid of frequencies are output 
by GALPROP either in Cartesian (l,b) or in HEALPix 
isopixelisation formafl The emissivity as seen by an 
observer at the solar position is also output as a function 

of (x,y,z,v). 

III. Radio surveys 

The aim is to describe both spectral and spatial prop- 
erties of the synchrotron sky, and not the fine angular 
details. For the spatial distribution the most useful data 
are the Haslam 408 MHz survey (|9|), since it has full 
sky coverage, a well-established calibration and zero 
level, and is dominated by synchrotron emission. The 
other data used in this paper are from 22 MHz: DRAO 
survey (02)); 150 MHz: flTT]; 1420 MHz: 02 and lfl3l . 
For WMAP frequencies we used the spectral-index maps 
generated from WMAP polarized data by |[T4l . scaling 
from 408 MHz to give total synchrotron intensity. 

IV. Procedure 

We have implemented many 3D models of the Galac- 
tic magnetic field in GALPROP in order to calculate the 
synchrotron emission from the Galaxy. Our approach is 
to use the most up to date models of the magnetic field 
from the literature (e.g. 03, 03, OS), 03, 03) and 
adjust the parameters in order to reproduce at the same 
time gamma-ray and radio observations. 

In future we will use a distribution of cosmic-ray 
sources and an electron spectrum that fits the Fermi 
gamma-ray spectrum and electron measurements. Then 
the value of total magnetic field will be adjusted to 
fit the synchrotron 408 MHz map. We will compare 
synchrotron latitude and longitude profiles with the 
available radio surveys. The best Galactic magnetic field 
model that fits simultaneously cosmic-ray, gamma and 
radio data will be presented at the conference. Results 
will be given using the cosmic -ray electron spectrum 
measured by Fermi [19|. 

V. Results 

Here we show an example of the outputs that can 
be obtained with GALPROP for a given model of the 
magnetic field. In this case the electron spectrum used 
was the so called conventional model defined in J6), 
based on pre-Fermi electron measurements. 

The regular Galactic magnetic field used in this ex- 
ample is an axisymmetric logarithmic spiral model with 
a local intensity of 6 and a pitch angle of 8.5°. 
The halo magnetic field is given by the z-dependence 
taken from [14|, but the parameters are adjusted to fit 
the latitude profiles of the synchrotron emission. The 
disk random field was taken as 3 /iG, constant in Galac- 
tocentric distance, with a halo component decreasing in 
Z with the same law as the regular one. Toroidal and 
dipolar magnetic field components such as in ifTTl are 
also added. 

'http://healpix.jpl.nasa.gov 



Synchrotron spectra for the inner and the northern 
Galaxy obtained with GALPROP are shown in Fig[T] 
compared with radio surveys. The spectrum is repro- 
duced except at the lowest frequencies where absorption 
is evident, which is not yet included in the model. Re- 
sulting latitude and longitude profiles of the synchrotron 
emission at 408 MHz are plotted in Fig. [2] compared 
with data from [9|. The general form of the profiles are 
already well fitted with this preliminary model. 




frequency, MHz 



Fig. 1. Synchrotron spectra for the inner Galaxy (left) (300°<1<60°, 
|b|<10°) and for the region 20°<1<180°, |b|<60° (right). Blue lines 
are obtained with the GALPROP code, while red points are radio 
surveys. The data are from 22 MHz: DRAO survey ([TOj); 150 MHz: 
QJJ; 408 MHz: 0; 1420 MHz: (TJ) and (TJ; WMAP: fl4l . 
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Fig. 2. Intensity profiles of synchrotron emission at 408 MHz in latitude for the inner Galaxy, |1| <60°, and the whole Galaxy, and in longitude 
for the Galactic plane region. Blue lines are obtained with the GALPROP code, while red lines are data from |9|. The zero level is corrected 
as described in [20] including subtraction of the CMB and extragalactic background. 



